use "Table2data.dta" 

*gen treat = 1 if total_intensity > .6
*replace treat = 2 if total_intensity > .5 & treat == .
*replace treat = 3 if total_intensity > .4 & treat == .

reg log_total_soilbase i.treat##i.year if year < 1935 
reg log_total_soilbase i.treat##i.year i.state if year < 1935 
